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Developing a theoretical framework to access quantum phases of itinerant bosons or fermions in two di- 
mensions (2D) that exhibit singular structure along surfaces in momentum space but have no quasi-particle 
description remains as a central challenge in the field of strongly correlated physics. In this paper we propose 
that distinctive signatures of such 2D strongly correlated phases will be manifest in quasi-one-dimensional 'W- 
leg ladder" systems. Characteristic of each parent 2D quantum liquid would be a precise pattern of ID gapless 
modes on the A^-leg ladder These signatures could be potentially exploited to approach the 2D phases from 
controlled numerical and analytical studies in quasi- ID. As a first step we explore itinerant boson models with 
a frustrating ring exchange interaction on the 2-leg ladder, searching for signatures of the recently proposed 
two-dimensional d-wave correlated Bose liquid (DBL) phase. A combination of exact diagonalization, density 
matrix renormalization group, variational Monte Carlo, and bosonization analysis of a quasi- ID gauge theory, 
all provide compelling evidence for the existence of a new strong-coupling phase of bosons on the 2-leg ladder 
which can be understood as a descendant of the two-dimensional DBL. We suggest several generalizations to 
quantum spin and electron Hamiltonians on ladders which could likewise reveal fingerprints of such 2D non- 
Fermi liquid phases. 



I. INTRODUCTION 

During the past decade it has become abundantly clear that 
there are several sub-classes of two-dimensional (2D) spin 
liquidsfi Perhaps best understood are the topological spin liq- 
uids that support gapped excitations which carry fractional 
quantum numbersi^^ii^ Another possibility are gapless spin 
liquids with no topological structure. These gapless spin 
liquids will generically exhibit spin correlations that decay 
as a power law in space and which can oscillate at par- 
ticular wavevectors. In such "algebraic" or "critical" spin 
liquid a'-^-^'^ these wavevectors are limited to a finite discrete 
set, and their effective field theories will often exhibit a rela- 
tivistic structure - as is the case when the spinons in a slave 
particle construction are massless Dirac fermions coupled to a 
U(l) gauge field. However, it is possible that the spin correla- 
tions exhibit singularities along surfaces in momentum space. 
These are analogous to the Fermi surfaces in a Fermi liq- 
uid but in such "quantum spin-metals" a quasiparticle picture 
will presumably be inapplicable. Spin liquids with a spinon 
Fermi surface are examples of particular interest, and have 
been studied in a number of RVB work a''^-'°'^'-'^-'-^ - most 
recently as a candidate for a spin liquid phas o^"^'^^'^^ in the or- 
ganic compound k-(ET)2Cu2(CN)3. 

Developing a theoretical framework for itinerant 2D non- 
Fermi liquids is arguably even more challenging than for spin 
liquids. Towards this end recent worki2, explored the possi- 
bility of an uncondensed quantum phase of itinerant bosons 
in two dimensions which is a conducting fluid but not su- 
perfluid. Because of the characteristic nodal structure ob- 
served by moving one boson around another, Ref. [13 called 
such phases "d-wave Correlated Bose Liquids" (DBL). Vari- 
ous physical correlators were found to be singular across sur- 
faces in momentum space in the DBL. For example, the boson 
occupation number n(q) is singular at some "Bose surfaces" 



like those illustrated in Fig. [T] Other examples with critical 
surfaces have been studied recently in Ref.USi 

The central idea underlying this paper is that 2D phases of 
quantum spins and itinerant fermions or bosons with singular 
surfaces in momentum space should have definite signatures 
when restricted to a quasi one-dimensional (ID) geometry, 
e.g. when the system is placed onto an A^-leg ladder. To be 
concrete, suppose we have a singular 2D surface in momen- 
tum space in the ground state of some square lattice quantum 
Hamiltonian. If we put the Hamiltonian on an A^-leg ladder, 
the transverse momentum qy becomes discrete. With large N 
one would expect that the short-range energetics that stabi- 
lizes the 2D phase will still be present on the ladder and will 
lead to a set of ID gapless modes where the discrete momenta 
cut through the singular surface, as illustrated in Fig.[T] 

The long wavelength physics of these ID modes will be 
described in terms of some multi-channel Luttinger liquid, 
behavior different than the parent 2D state. But the number 
of gapless ID modes and their respective momenta character- 
izing the oscillatory power law decays would provide a dis- 
tinctive "fingerprint" of the parent 2D liquid. For example, 
the number of gapless ID modes would grow linearly with 
N and the sum of their momenta could satisfy some gener- 
alized Luttinger volume. Moreover, these signatures could 
be potentially exploited to approach the 2D phases from con- 
trolled numerical and analytical studies in quasi-lD. Indeed, 
as we demonstrate explicitly in this paper, the (putative) sin- 
gular surface of the 2D DBL phase is already manifest on the 
2-leg ladder! 

The itinerant uncondensed DBL of Ref. [l^ is constructed 
by writing a boson in terms of fermionic partons, 

= 44 : (1) 

and considering un-coupled mean field states of the di and 
^2- The DBL phase on a square lattice is obtained by tak- 
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FIG. 1: Schematics of discrete qy — 2n/N cuts throughi Bose sur- 
faces in some anisotropic 2D system placed on an A^-leg ladder; here 
for illustration N — 8. Such Bose surfaces appear in a DBL con- 
struction described in the text; to avoid clutter, we do not show the 
Fermi surfaces of the di and d2 slave particles. 

ing the di fermions to hop preferentially along the x-axis and 
the d2 fermions to hop preferentially along the y-axis. One 
can enforce the square lattice symmetry in the boson state by 
requiring that the di and d2 hopping patterns are related by 
a 90 degree rotation. To recover the physical Hilbert space 
requires the condition d\d-^ = d\d2 = b'^b at each site. 
Within a Gutzwiller projection approach the corresponding 
boson wavefunction is taken as product of two fermion Slater 
determinants, 

'^hos{ri,r2,...) = *di(r-i,r2,.. .) • ^^^(ri, ra, . . . ) . (2) 

In the DBL the di and ^2 Fermi surfaces are different, com- 
pressed in the x and y-directions, respectively. This results 
in a characteristic d-wave-like nodal structure observed when 
moving one boson around another in 4'bos- The boson mo- 
mentum distribution function n(q) in the DBL is singular on 
two surfaces that are constructed from the di and ^2 Fermi 
surfaces as enveloping kpi ± kp2 surfaces. More details are 
in Ref. 17, while Fig. [T] shows one example of the kpi± kF2 
loci for an open di Fermi surface and a closed d2 Fermi sur- 
face (the Fermi surfaces themselves are not shown; this DBL 
state would be anisotropic on the square lattice but is closer 
to the ladder states considered here where there is no 90 de- 
gree rotation symmetry). As an alternative to the Gutzwiller 
construction one can project into the physical Hilbert space 
within a gauge theory approach, coupling di and d2 with op- 
posite gauge charges to an emergent U(l) gauge field. 

Ref. '17' also proposed a frustrated boson model with ring 
exchanges that can potentially realize such DBL phases: 

H = Hj + idling , (3) 

= E (fei^r+A + H-c), (4) 

r; (1—x.y 

iJring = KY^{blb,+A+,+yK+y + ^-C-) , (5) 



with J,K > Q. In addition to the usual boson hopping term 
this Hamiltonian contains a four-site ring term acting on each 
square plaquette. With K positive this term is "frustrating", 
violating the Marshall sign rule of the hopping term - mak- 
ing the system intractable by quantum Monte Carlo simula- 
tions. This boson Hamiltonian was constructedii by taking 
the strong coupling limit of the lattice gauge theory for the 
di and d2 fermions. Increasing the disparity in the fermion 
hopping along the x and y axes corresponds to an increase in 
K. 

In this paper, we make a first step in the proposed program 
of ladder studies by exploring the ring model Eq. (O on a 2- 
leg ladder. Of course, the picture of discrete qy cuts through 
two-dimensional surfaces is pushed to the extreme here. Nev- 
ertheless and quite remarkably, in a wide parameter regime of 
intermediate to large ring couplings, we find an unusual phase 
which can be understood via the DBL construction Eq. (|2|i. 
This is a strong-coupling phase with no perturbative picture 
in terms of the original bosons. The slave-particle approach 
provides a new starting point, and the resulting gauge theory 
of the DBL can be solved by the conventional ID bosoniza- 
tion toolsr2i22i2iiS providing a consistent picture of the new 
phase. 

The paper is organized as follows. Section sets the stage 
by describing DBL states on the 2-leg ladder The main Sec- 
tion|lIl]presents results of exact numerical studies of the J—K 
model. The focus here is on the DBL phase that emerges 
when the conventional "q = liquid" is destroyed by the ring 
exchanges. However, the phase diagram is even more rich 
and contains two more unusual phases in which bosons are 
paired; one of these phases can be accessed analytically as an 
instability of yet a different DBL phase. Appendix|A] summa- 
rizes the technical gauge theory description of the DBL phases 
while Appendix IB] presents simple analytical results for the 
ring model in some limiting cases. SectionHVlconcludes with 
an outline of future program of ladder studies. 



II. DBL STATES FOR HARD-CORE BOSONS ON A 2-LEG 
LADDER 

The DBL construction Eq. proceeds by taking distinct 
hopping ground states for the di and d2 fermions. The hop- 
ping patterns of the two fermion species can be independent 
and need to only respect the ladder symmetries. There is a 
bonding and an antibonding band for each species. We take 
the di fermions to hop more strongly along the chains and al- 
ways assume that the chemical potential crosses both bands 
as illustrated in Fig.|2] We take the d2 fermions to hop more 
strongly on the rungs and consider two cases shown in the bot- 
tom panels in Fig. |2] The "DBL[2,1] phase" when only the 
bonding band of the d2 fermions crosses the Fermi level, and 
the "DBL[2,2] phase" when both ^2 bands are partially oc- 
cupied - the first and second argument in the square brackets 
refer to the number of involved di and d2 bands respectively. 

Let us first discuss properties of the DBL[2,1] phase. We 
denote the Fermi momenta of the di bonding and antibond- 
ing bands as k!'p\ and fc^j^, corresponding to the transverse 
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FIG. 2: (color online) di and d2 bands for the DBL construction 
on the 2-leg ladder. Top panel: di fermions partially occupy both 
the bonding and anti-bonding bands. Bottom left: In the DBL[2,1] 
phase, d2 fermions occupy only the bonding band. Bottom right: In 
the DBL[2,2] phase, d^ fermions partially occupy both bands; this 
state is likely unstable towards a boson-paired phase, as discussed 
in Sec. lA 31 In the wavefunction construction, only the occupied 
orbitals and not the band dispersions are important. 



momenta ky = and ky = tt respectively (which is a con- 
venient alternative to the quantum number for the leg inter- 
change symmetry). We denote the Fermi momentum of the ^2 
bonding band as fcF2 (there is no d2 antibonding occupation 
in the DBL[2,1]). The Fermi momenta satisfy kpl + k'^pl = 
kp2 ^ "2,71 p, where p is the boson density per site and remem- 
bering that there are two sites in each rung while the lattice 
spacing is set to one. 

In a mean field treatment, the boson Green's function is a 
product of two fermion Green's functions and therefore os- 
cillates at different possible kpi ± ^'F2 wavevectors and de- 
cays as l/x^. Going beyond the mean field as described in 
Appendix |A1 we expect that the oscillations at wavevectors 
q^^^^ = kp2 — kpl^'' have slower power law decay while 

the other two wavevectors kF2 + kp(^^ have faster power law 
decay than the mean field. The complete result is, dropping 
all amplitudes. 
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The first line shows the enhanced correlators. These can be 
guessed using an "Amperean rule" mnemonic inherited from 
the fermion-gauge studies in 2D and verified by the bosoniza- 
tion solution in Appendix|A| Composites that involve fermion 
bi-linears with opposite group velocities but that produce par- 
allel gauge currents are enhanced by the gauge fluctuations. 
Thus, low-energy "bosons" df/'^^^ dlj^ carrying q^^^^ are ex- 
pected to be enhanced. The contributions in the second line, 
on the other hand, can be shown to be always suppressed be- 
yond the mean field. Appendix lAl presents a complete low- 



energy theory of the DBL[2,1] which - once the gauge fluctu- 
ations are treated - contains two phonon modes and in princi- 
ple allows calculating all exponents. Since the power laws are 
affected also by short-range interactions, and there are many 
parameters allowed in the theory, we do not pursue this calcu- 
lation in general, but mainly rely on the Amperean rules that 
capture the gauge fluctuation effects; this is reasonable if the 
gauge interactions dominate. 

The boson singularities at non-zero wavevectors become 
manifest upon Fourier transformation of Gb (r) to obtain the 
boson occupation number n(q). For example, at the enhanced 

wavevectors {qi?'' , 0) and (gi'^' , tt), we have 



n{qi''y'^ + Sq^,ky) - \Sq^\^ 



(8) 



The boson density correlations are analyzed similarly. In a 
naive mean field, -D&(r) = (n(r)n(O)) ~ Ddi{r) + Dd2(r), 
which oscillates at various "2k vectors and decays as l/x^. 
Beyond the mean field, the power laws get modified: 
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The distinct momenta are: (2fc^°^,0), {2k''p^,0), (2fcF2,0), 

{kfl + k'^p^TT), {k^pl - k'p^,7T), and (0,0). The first four 
involve a particle and a hole of the same species but with op- 
posite group velocities and are therefore expected to be en- 
hanced by the gauge fluctuations; in the above equations, this 
is indicated schematically with l/a;^^' . However, these cor- 
relators can also be affected - either enhanced or suppressed 
- by short-range interactions, and without knowing all param- 
eters in the theory, we cannot calculate the exponents reliably. 
The oscillation at wavevector (kpl — k'^p^ , tt) can be shown to 
be always suppressed beyond the mean field. Finally, the zero 
momentum power law remains unchanged. 

More details on the DBL[2,1] -phase can be found in Ap- 
pendix [A] In particular, we show that the two-boson Green's 
function exhibits some internal d-wave character, which origi- 
nates from the non-trivial wavefunction signs. It is these signs 
that prompted the name "D-wave" Boson Liquid (DBL) in the 
2D continuum setting,'^ where the wavefunction goes through 

a sequence of signs H h— upon taking one particle around 

another. However, such correlations do not necessarily mean 
that the system is near a d-wave-paired phase. For example, 
we also examine a potential instability of the DBL[2,1] driven 
by an allowed non-marginal four-fermion term and find that 
the resulting phase is a boson-paired liquid with an internal 
s-wave character 

Let us briefly mention the DBL[2,2] case, where in the 
mean field both di and d2 have partial occupations of both 
the bonding and antibonding bands. One can calculate vari- 
ous correlations in the mean field as before, and the structure 
is more rich since there is an additional band present. Some 
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FIG. 3: (color online) Phase diagram for the 2-leg boson system 
at filling p = 1/3. The symbols represent DMRG estimates with 
shaded grey circles indicating the superfluid phase, the filled blue 
circles the d-wave correlated Bose liquid (DBL), and the open red 
diamonds the s-wave paired state. The schematic line boundaries are 
obtained from VMC calculations. 



FIG. 4: (color online) Phase diagram for the 2-leg boson system at 
filling p = 1/9. The symbols indicate DMRG estimates for the var- 
ious phases with shaded grey circles for the superfluid phase, filled 
blue circles for the DBL, and open green diamonds for the d-wave 
paired state. The schematic line boundaries are obtained from VMC 
calculations. 



details are given in Appendix I A 3! Beyond the mean field but 
ignoring non-marginal four-fermion terms, this phase would 
have three gapless modes. However, an analysis of allowed 
interactions suggests that this phase has a strong instability 
towards a boson-paired phase with an s-wave character. Re- 
markably, as we will see in the next section, a boson-paired 
phase with similar properties is found in the ring model in 
the regime of small interchain coupling Jj^ where we initially 
hoped to find the DBL[2,2] phase. Thus, the slave-particle 
formulation and the gauge theory description solved via sub- 
sequent bosonization open a non-perturbative access to this 
interesting boson-paired phase. 



III. FRUSTRATED J-K MODEL ON THE 2-LEG LADDER 

We study the ring model Eq. ([3j on a 2-leg ladder with a 
boson hopping term J that allows different amplitudes along 
and perpendicular to the chains: 

= -jY.if'lK+, + H.c.) - ^± EK.,i)^..2) + H-cXl 

r X 

where r = {x, y) represent integer lattice sites. We take y = 
1, . . . ,Ly with Ly ^ 2 and x ~ 1, . . . ,Lx with the length 
of the 2-leg system. The ring terms are associated with the 
square plaquettes of the 2-leg ladder 

We use exact diagonalization (ED) and density matrix 
renormalization group (DMRG)2^i2i methods supplemented 
with trial wavefunction variational Monte Carlo (VMC}'^^ to 
determine the nature of the ground state of the Hamiltonian 
Eq. ([3]) at boson filling number p = jT"' where Nb is the 
total number of bosons in the system. The obtained phase di- 
agrams at densities p = 1/3 and 1/9 are shown in Fig. [3] and 



Fig.|4]respectively. These are typical results for the boson ring 
model away from half-filling. 

Using ED and DMRG, we find four different quantum 
phases upon increasing the ring-exchange coupling K and 
varying the interchain coupling J±. At small K, the "super- 
fluid" phase (SF) with q = quasi-long-range order (QLRO) 
is stable for both filling numbers. The DBL phase, which is 
produced by the ring exchanges, dominates the intermediate 
parameter space at higher filling number p ~ 1/3, while it oc- 
cupies a substantially smaller region at p = 1 /9. Interestingly, 
at p = 1/3 and small J±, an s-wave paired phase emerges 
between the SF and DBL as a consequence of the competi- 
tion between the boson hopping and ring-exchange terms. At 
p = 1/9 and small J±, a d-wave paired phase is found ad- 
jacent to the DBL. Finally, for strong ring exchanges phase 
separation eventually wins near iiT ~ 4 J for both fillings, sep- 
arating into a region with p = 1/2 and an empty region p — 0. 
The characteristic features of each phase will be discussed in 
the following sections. 

Our complementary VMC study is centered around the 
DBL phase. The DBL wavefunction Eq. Q is constructed 
from the slave-particle formulation and thus allows us to com- 
pare between the exact ED/DMRG results and the gauge the- 
ory description. Jastrow SF and boson-paired trial wavefunc- 
tions (described in Appendix |C]l are also considered in order 
to better understand the exact numerical results. We note 
that this variational approach has been rather successful in 
one dimension, e.g., in the original studies of the ID i — J 
model] We borrow some ideas from these studies such 

as allowing the fermionic determinants to also play a Jastrow 
role via dct sign(det)|det|^, where p is a variational pa- 
rameter. In our ring model, the phase boundary estimated 
from VMC almost coincides with the one determined from 
DMRG as indicated in Fig. [3] for p = 1/3 filling. At smaller 
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filling number, DMRG and VMC also find the same topology 
of the phase diagram as given in Fig.H] 

We first describe measurements in the ground state. Im- 
posing periodic boundary conditions in the i-direction in our 
numerical calculations, we can characterize the ED states by 
a momentum fc^. In our DMRG calculations we work with 
real-valued wavefunctions which gives no ambiguity when 
the system exhibits a unique ground state that caries zero mo- 
mentum. If the ground state carries nontrivial momentum k^, 
its time-reversed partner carries —k^, and the DMRG state 
is some real-valued combination of these. While the lattice- 
space measurements may depend slightly on the specific com- 
bination (but vanishingly in the large system limit), the mo- 
mentum space measurements n(q) and D{q) described below 
are unique and insensitive to this. In our DMRG calculations, 
more than 1500 states were kept in each blocks^ to ensure 
accurate results, and the density matrix truncation error is of 
the order of 10"^. The typical error for the ground-state en- 
ergy is of the order of 10^^ for all systems we have studied. 
The relative error in the correlation functions varies between 
10~^ — 10~^, depending on the type of the correlations and 
the spatial distance between operators. 

We measure boson correlator in the ground state 



Gfc(r,r') = (6t6r,) 
and its Fourier transform 
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FIG. 5: (a) The boson occupation number 7i(q) and (b) the density- 
density structure factor -D(q) for the 2-leg system at p = 1/3, 
Jx = J, K — J; the system length is Lx ~ 48. The results are 
representative for the superfluid phase. 



which is interpreted as mode q occupation number 
We also measure the boson density correlator 

A(r,r') = ((pfc(r)-p)(p,(r')-p)), 



(14) 



where p is the average boson density. The density structure 
factor is 



A(q) 



1 



LxLy , 



A(r,r')e 



iq-(r-r') 



(15) 

For both the boson or boson density correlations, a power law 
1/xP envelope in real space gives rise to a singularity 
in momentum space. 

Finally, we measure two-boson (pairing) correlator, 

P26(ri,r2;rl,r^) = {bHri)b\r2)b{r[)b{r'^)) , (16) 

which creates a pair of bosons at ri, r2 and removes a pair at 
r[ , r'^- This is useful for detecting boson-paired phases where 
the single boson correlator decays exponentially. In the pair- 
ing correlator figures below, we present P26 for a fixed cre- 
ated +45° diagonal pair {ri = (0,1), V2 = (lj2)} near 
the origin, while the removed pair is either a +45° diago- 
nal {v\ = (x, 1), v'2 = {x + 1,2)} or a -45° diagonal 
pair {r'j = (a;, 2), r'j = (x + 1,1)}. We have measured 
the correlator for other pair orientations as well, but find that 
the diagonal-diagonal correlations are the most distinguishing 
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FIG. 6: The boson pairing correlations in real space for the SF sys- 
tem described in Fig. [5] We show correlations for pairs on the di- 
agonals d and d, but roughly similar values are obtained for all pair 
orientations. 



ones in the two paired phases. When the amplitudes for the 
±45° diagonals have the same sign, we refer to this as s-wave 
character; on the other hand, when the signs are opposite, we 
call this d-wave. An extended discussion can be found in Ap- 
pendixlAlfollowing Eq. ( IA13l l. 
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A. Superfluid Phase 

Let us first discuss the phase at very small K. The system is 
not frustrated in the K = limit and the ground-state wave- 
function (WF) components in the boson occupation number 
basis should all be of the same sign (positive) to gain the 
lowest hopping energy. Interestingly, the ED calculations for 
small Lx ~ 12 and L^; = 18 systems find that this feature of 
the WF remains robust up to a finite value of K ^ 2 J, close 
to the stability boundary of the SF phase. This is illustrated in 
the medium panel of Fig. [7] where the average sign of coeffi- 
cients in the ground-state wavefunction is plotted. The boson 
occupation number n{qx , Qy) has a sharp peak at the wavevec- 
tor q ~ (0, 0). This is shown in Fig. |5](a) for a system at 
p = 1/3, J± = J, K ^ J, and length = 48, calculated us- 
ing the DMRG method. For this system size, the peak value at 
q = (0, 0) is about 100 times larger than the value at wavevec- 
tor q = (0, tt). This is consistent with the expected singular 
behavior n{qx,0) ~ k^l''''^^ and non-singular n{qx,Tr) in 
the SF phase with ID QLRO; the boson correlator decays as 
l/x^'' in real space and we estimate pi, w 0.25. 

As shown in Fig. |5] (b) for the same system, the boson 
density-density structure factor D{qx,0) has a \qx\ depen- 
dence around qx = reflecting the global boson number con- 
servation, while D{qx,7T) is nonzero and has no visible singu- 
larities. 

The pairing correlations in real space show power-law be- 
havior, with nearly equal magnitudes for all pair orientations. 
This is shown in Fig. |6] where we create a +45° diagonal 
pair near origin and remove either a +45° diagonal pair or 
a —45° diagonal pair near x. Furthermore, the exponent p2b 
of the power-law 1 / x^^'' in the pairing correlator is found to 
be p2b ~ 1, which is about four times that for the single boson 
correlator, p2b ~ 4pf,. Such behavior of the pairing corre- 
lations is expected in the SF phase, which has QLRO in the 
single boson correlations. 

These features characterize the SF phase. It also has a finite 
PsL, where ps is the "superfluid stiffness", which we can mea- 
sure in our numerical calculations by imposing a small twist in 
the boundary conditions. Results from ED of a small system 
with Lx — 12 are given in the top panel of Fig.|2]and related 
DMRG measurements for larger systems are in Fig.[TTl these 
are discussed in more detail in Sec. lIIIBl 

B. D-Wave Correlated Bose Liquid (DBL) 

As we further increase the strength of the ring-exchange K, 
the system undergoes a (first order) phase transition where the 
stiffness ps drops steeply by a factor of two, see Figs. 1711 II 
and the subsection below. There is also an abrupt change in 
the character of the wavefunction sign which we can monitor 
in the ED analysis for small system sizes by simple counting 
and shown in the middle panel in Fig.|2l^^ This new phase has 
unusual phenomenology and can be identified as the DBL[2, 1 ] 
phase: There is no boson peak at q = (0, 0), but instead the 
boson occupation number n{qx, qy) shows singular points at 

nonzero wavevectors (igi"' , 0) and (igi'^' , tt), as if there is 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 




ring exchange K / J 



FIG. 7; Exact diagonalization results for a system of size Lx = 12 at 
filling p = 1/3 and J± — J. Shown are the superfluid stiffness (top 
panel), the average sign of the ground-state wavefunction (medium 
panel) and the average double occupancy of rungs (lower panel) as a 
function of the ring exchange K. 



a "condensate" of bonding and antibonding boson modes at 
these qx- Fig. [8]presents characteristic plots of the n{qx, qy) 
and the density-density structure factor D{qx, qy) deep in this 
phase for filling p ^ 1/3 taken at J± = J, K = 3 J, and a 
system length Lx ~ 48. 

The singular boson wavevectors q^^ and qi^'^ in general 
vary with parameters K and J±. An example is shown in 
Fig. |9]for the system at p = 1/3 with J± = J and vary- 
ing K. Just into the DBL phase near L\ > 2, we find a 
small q^x'^ and a large qi^K In the wavefunction construc- 
tion, this corresponds to rather different k^pl and k^pi, with 
significantly more di fermions in the bonding than antibond- 
ing bands, see Fig. [J] When we increase the ring coupling, 
decreases while gi'^' increases, corresponding to moving 
the di fermions from the bonding to antibonding levels and 
making the two bands more equally populated. Interestingly, 
we find that the relation gi"' + qi^'^ = 2Trp is satisfied, consis- 
tent with the prediction from gauge theory of a Luttinger-type 
theorem in the correlated boson phase. 

The identification with the DBL is bolstered by a trial wave- 
function study. The optimal such state for the system in Fig. [8] 
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FIG. 8: (a) The boson occupation number n(q) and (b) the density- 
density structure factor D{q) for the 2-leg system at p — 1/3, 
J± = J, K = 3J; the system length is Lx = 48. The results 
are representative for the DBL phase. 
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The length of the system is = 48 and the total boson num- 
ber is Nb = 32. The di determinant is constructed from 
n[^^ = 20 bonding and n[^^ = 12 antibonding orbitals, 
while the d-z determinant is constructed from bonding orbitals 
only. To define the orbitals, we use antiperiodic boundary con- 
ditions ("abc") along the x-direction for both di and ^2, which 
gives closed shells for the above occupation numbers, while 
the physical boson wavefunction respects periodic boundary 
conditions and carries zero total momentum. The determinant 
powers are understood as (det)^' = sign(det)|det|''. The d2 
determinant can be written explicitly as 



(0) _ 



= 0] oc JJ^ sm 



Tr{xi — Xj) 



(19) 



and in particular prevents two bosons from sharing the same 
rung. In fact, the exact ground states in the ED/DMRG have 
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FIG. 9: Positions of the singular momenta in the DBL phase as a 
function of the ring exchange K for a system with J± = J and 
filling p = 1/3. Data for systems with length Lx = 36 and Lx = 48 
are given. The momenta are determined from the peak values of 
n{qx, Qy) such as shown in the top panel of Fig.[8] 



very small double occupancy of rungs deep in this phase, as il- 
lustrated in the lower panel of Fig.|2l The characteristic boson 



(0) 



wavevectors in the above DBL state are Qx = kf2 

127r/48 and gi''^ = A:f2 - k^pi = 207r/48. These bosons 
are created by occupying di and d2 orbitals on the opposite 
sides, which is motivated by the Amperean enhancement of 
such composites as described earlier. The overall match of the 
trial wavefunction results with the exact DMRG correlations 
is striking, reproducing the singular features as illustrated in 
Fig.E 

Turning to the density structure factor D{q.j.,qy), there are 
many features but they are also weaker than in the n{qx, qy). 
For qy = 0, we expect a linear | g^; | -behavior near qx = 

and signatures at 2fc^"^' = 407r/48, 2k^p^ = 247r/48, which 
we see, and also at |2fci?2| = 327r/48 mod 27r, which we 
do not see in either DMRG or VMC (we can understand this 
in the wavefunction because the power for the det2 is small 
and suppresses density correlations coming from this piece). 



At ( 
and at k 



77, we expect signatures at kpl + k 



(0) 
Fl 



Fl 



3277/48 

877/48, both of which we see. As dis- 



cussed in Appendix|Al it is a complicated matter to predict the 
exponents of the different singularities, and given our limited 
system size we do not attempt measuring these. It is remark- 
able that the DBL wavefunction can reproduce the singular 
features, but we do not know at this point whether it can re- 
produce all exponents in general. 

Finally, in Fig.[TO]we show the diagonal-diagonal pairing 
correlator. It is more short-ranged than in the SF phase as well 
as the two paired phases described below, but is characteristic 
of the DBL phase. First, P2b{x) has robust sign oscillations at 
wavevector 277/?. Second, the +45° and —45° diagonal pairs 
have opposite signs; the nontrivial sign structure in the wave- 
function gives some "D-wave"-correIated character to the bo- 
son liquid, which is observed in the two-boson correlator. 
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FIG. 10: The boson pairing correlations in real space for the DBL 
system described in Fig.[8]calculated by DMRG and VMC. We show 
correlations for pairs on the diagonals d and d, while values for other 
pair orientations are much smaller. 



Detection of Boson Bonding and Antibonding Propagation 
by Two Stiffnesses 

In this subsection we want to further highlight some char- 
acteristic difference between the SF and DBL phases. In the 
SF state, the wavefunction has QLRO at q = 0. Imposing 
a twist boundary phase 9 in the boson hopping term in the 
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Hamiltonian (|3]l, one can measure the superfluid stiffness as 
Ps = d^E/dd^, which is found to scale as l/L^ in the SF 
phase. This stiffness calculated from ED is shown for the 
small Lx = 12 system in the top panel of Fig.|2l For larger 
system lengths, the stiffness can alternatively be estimated by 
DMRG calculations of the energy difference for a9 = tt twist, 
i.e., Ps - AS = E{e tt) - E{e = 0). 

In the top panel of Fig. [TT] we show L^^E as a func- 
tion of the ring exchange K/J for several system sizes = 
24, 36, 48. For isotropic hopping J± = J, we find that L^AE 
is essentially constant in the SF phase, but drops abruptly at 
the SF to DBL transition. In the DBL phase, L^AE shows 
rather irregular behavior including sign changes when vary- 
ing either K or L^- For a fixed system length L^, the more 
abrupt changes can be matched with the step-like changes 
of the singular momenta positions shown, e.g., in Fig. |9]for 
= 36, 48. For fixed ring exchange K, the abrupt changes 
with varying are indicative of strong incommensurate cor- 
relations, where different offer varying degree of matching 
at the (twisted) periodic boundary: AE > means that 9 — 
boundary conditions provide better matching than 9 ~ tt, 
while the situation is reversed if AE < 0. 

A finite |Lj;Aii^| in the large system limit in the above set- 
ting implies that there is a propagating "bonding" mode (i.e., 
with transverse momentum qy = 0). Based on our picture of 
the DBL phase, there must also be a propagating "antibond- 
ing" mode with qy — tt, which we would like to detect directly 
and contrast with the SF phase where there is no such mode. 

To this end, we design a system composed of two halves 
connected at both ends, where the left half is a {J, J± = 
J, isr)-HamiItonian while the right half is a (J, J± ~ — J, K)- 
Hamiltonian. First, we note that a model with J± ~ —J 
can be related to a model with J± = J by a transformation 
b{x,y) { — l)yb{x,y). This ensures that the boson density 
is the same in the two parts. Crucially, this transformation 
interchanges the notions of bonding and antibonding modes. 

If the designed system is in the regime of the regular SF 
phase, the gapless bonding bosons in the half with J± — J 
cannot penetrate into the other half where the gapless mode is 
anti-bonding, and vice versa. Since the bosons are not able to 
propagate around the loop, a 9 = n twist should not produce 
significant change AE which will vanish exponentially with 
Lx- On the other hand, in the DBL phase there are gapless 
modes of both bonding and antibonding type in either region, 
so that bosons can propagate around the loop, and the 7r-twist 
should change the energy by an amount proportional to l/L^- 
The "bonding/antibonding propagation stiffness" L^AE for 
the designed two-half system is shown in Fig.fTTTb) as a func- 
tion of K/J. Indeed, we find a vanishingly small L^AE 
in the SF phase and a rather irregular but nevertheless finite 
L^AE in the DBL phase.-^ 



FIG. 1 1 : Detection of the boson propagation around the loop by 
applying a phase twist 9 to the periodic boundary connection along 
the f-direction. The energy difference AE — E(9 = tt) — E{9 = 
0), multiplied by Lx, (top panel) for a uniform system with Jx = 
J for all bonds; (lower panel) for a specially designed system with 
Jx = J in one half and Jx = — J in the other half of the system. 



C. Extended S-Wave Paired Boson Phase 

For our 2-leg ladder system one might expect to find the 
DBL[2,2] phase if the interchain coupling J±^ is much smaller 
than the coupling J along the chains: When the boson hop- 
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momentum q 

FIG. 12: (a) The boson occupation number n(q) and (b) the density- 
density structure factor -D(q) for a 2-leg system at p = 1/3, J±_ = 
O.IJ, K — 1.4J; the system length is Lx = 48. The results are 
representative for the s-wave paired phase. 



The observed phase can be identified as an s-wave bo- 
son paired phase with QLRO pairing correlation as shown 
in Fig. [13] The pairing correlation P2b{x) in this phase is 
strongly enhanced at large distances x when compared to the 
correlation in the DBL phase; the values are roughly similar 
for pairs straddling the chains but smaller for pairs inside the 
chains. The decay of P2b{x) can be fitted by a power law 
behavior-- 1/\x\P''>' + l/\L^ - x^""" withp2b « 1- 

We have also checked the single boson correlation in the 
real space: It is much reduced compared to the nearby SF and 
DBL phases and decays very quickly - essentially exponen- 
tially; in particular, it is much smaller than the boson-pair cor- 
relation beyond few lattice spacings. We then conclude that 
there is a gap to single-boson excitations, but not to pairs. This 
identification of the phase is also supported by a VMC study 
using paired boson wavefunctions described in Appendix IC 21 
which improve the energetics over that of the trial DBL states 
and reproduce qualitatively the DMRG correlations. 

Interestingly, the s-wave paired phase can be regarded as 
an instability of the DBL[2,2] phase. In Appendix lA 31 we 
argue that without specially adjusted short-range interactions, 
the DBL[2,2] is unstable once the gauge fluctuations are in- 
cluded, and the most natural outcome is a paired phase with 
the same phenomenology as described above. In particu- 
lar, the pairs are s-wave, carry zero momentum, and show 
QLRO with power law x"'^^'', while we also predict power 
law charge correlations at wavevector {2Tip, tt) with x^'Pp . In 
Appendix lA 31 we also predict that the two exponents are in- 
verses of each other, P2bPp — 1- The DMRG estimates of 
the power laws differ somewhat from this, but they are hard 
to make accurate given the slow decay of the coiTelations and 
our limited system sizes. 



ping between the chains is reduced, one expects that the 
fermion hopping between the chains is also reduced for both 
species, which in turn could bring the antibonding band of the 
fermions to cross the chemical potential as shown in the 
bottom right panel of Fig. |2] DMRG studies of the ring model 
at filling p = 1/3 and small < 0.3 indeed find a phase 
that is distinct from both the DBL[2,1] and SF as shown in 
the phase diagram Fig. [3] around K ^ lAJ. However, our 
measurements indicate that this phase is not a DBL[2,2], but 
a paired phase. As shown in Fig.[T2]for a system at p = 1/3 
with = O.IJ, K = 1.4 J, and length = 48, this phase 
is characterized by a quite flat region in the boson occupation 
number n{qx,qy) near = without any singular points, 
which is different from both the SF and DBL phases. On the 
other hand, the density-density structure factor D{qx, Qy) has 
a singular behavior at a wavevector {2ttp,tt). This singular 
behavior can be approximated by ~ with pp be- 

tween 1.33 — 1.50 providing the best fit to the data, which 
coiTesponds to power law envelope x^^'' in real space. [Note 
that when J± ^ the number of bosons in each chain is in- 
dependently conserved and D{qx, it) must vanish at q^ = 0, 
which explains the small value of D{qx = 0, tt) extracted by 
DMRG when Jj, = O.IJ in Fig. [12]] 



D. D-Wave Boson Paired Phase 

Consider now bosons at small filling p = 1/9, Fig. [4] As we 
further increase K, the DBL phase in the intermediate range 
becomes unstable and a new phase emerges. As shown in 
Fig. [14] for a system at filling p = 1/9 with J± = O.IJ, 
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FIG. 13: The boson pairing correlations in real space for the same s- 
wave paired system described in Fig.[T2] We show diagonal pairs but 
roughly similar correlations are obtained also for pairs on the rungs, 
while correlations for pairs in the chains are significantly smaller. 
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FIG. 14: (a) The boson occupation number n(q) and (b) the density- 
density structure factor -D(q) for a 2-leg system at p = 1/9, J± = 
O.IJ, K = 3 J; the system length is ~ 54. The results are 
representative for the d-wave paired phase. 



K = 3 J, and length = 54, this phase is characterized by 
n{qx, 0) ~ n{qx, tt), which is satisfied very accurately and is 
indicative of strongly suppressed boson correlator between the 
chains. We think that the peaks in n{q) are non-singular and 
represent short-range boson correlations, while G'f,(r) decays 
exponentially on long distances: While we cannot ascertain 
the exponential decay with our system sizes, the single boson 
correlator is already more than an order of magnitude smaller 
and decays much faster than the pair correlator described be- 
low. This phase is further identified as a d-wave boson paired 
phase with power-law pairing correlation as shown in Fig.fTSl 
The diagonal pairing correlations P2b{x) show opposite signs 
depending on the relative orientation of the two pairs: positive 
when both are +45° diagonals and negative when one is +45° 
and the other —45°, which is consistent with the d-wave pair- 
ing picture. We find power law correlation P2b{x) ~ l/x^^'' 
with exponent p2b ~ 0.85 

The identification of this as d-wave paired phase is sup- 
ported by the reasoning at low densities described in Ap- 
pendix lB3] Thus, for a pair of bosons on an otherwise empty 
ladder, there is a transition at K/ J = 2 where it becomes fa- 
vorable for bosons to form a molecule with internal d-wave 
character. If we now have a system of bosons at small den- 
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FIG. 15: The boson pairing correlations in real space for the same 
system described in Fig.ll4l 



sity, it is natural that there will be a parameter range where 
the bosons are paired and the resulting molecules form a Lut- 
tinger liquid. 

To bolster this identification, we also studied such states 
variationally. As described in Appendix IC 21 we construct a 
trial wavefunction as a product of a Pfaffian and a determinant, 
since these are simple to work with in VMC. Optimal trial 
state for the system in Fig.[T4lis 



*1: 



= Pfi X 



12, N^'^^ = 0] 



0.3 



(20) 



The system length is ~ 54 and the boson number is Ni, = 
12. The Pfaffian is for an antisymmetric (spinless fermion) 
pair-function 5(2:,?/; x',?/') = (5y^y'Sign(j/ - j/')e~l^~^''l/? 
that connects only sites on different chains; the optimal "size 
of the pair" is ^ = 1.8. This is a "py" pairing and when com- 
posed with a "pa:" character present in the determinant this 
gives a d-wave character to the boson pairs. Strictly speaking, 
the above state is valid only in the = limit, since it pre- 
vents fluctuation of the boson number in each chain. Still, as 
we can see from Fig. [141 the variational wavefunction repro- 
duces the DMRG correlations surprisingly well in this phase. 
The boson momentum distribution function for the trial state 
is (by construction) non-singular, so that the close agreement 
with the DMRG found for Lx = 54 strengthens the identifi- 
cation of this phase as a d-wave paired. 

Finally we mention the signatures in the boson density 
shown in the lower panel in Fig. [141 D{Qx, 0) shows the fa- 
miliar \qx\ at small momenta and also a singularity at 27rp. 
The latter gives a wavelength equal to the spacing between 
the pairs, which is expected in a Luttinger liquid of such 
molecules. The density singularity is not too strong and is 
in line with the slow decay of the pair correlator. 



E. Phase Separation 

As we further increase the ring exchange towards K > 4 J, 
the uniform phase becomes unstable to phase separation. For 
such large K, the boson density will separate into half-filled 
and empty regions, which is energetically favored by the ring 
exchange term. In the ED of the small La, = 12 system. 
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phase separation is identified to happen at K > 5J where all 
the lowest energy states from different momentum sectors be- 
come nearly degenerate. In the DMRG calculation, the phase 
separation is indicated by a nonuniform boson density in the 
obtained ground state, which breaks the translational symme- 
try due to small perturbations from the cut-off of the Hilbert 
space in DMRG. The phase separation occurs around K > 4 J 
as observed from DMRG. Finally, the phase separation can 
also be studied in VMC by measuring trial energies for dif- 
ferent boson densities and performing Maxwell construction 
(this is described in more detail, e.g., in Ref. J_2 Sec. Vll for 
the 2D ring model). 

The tendency towards phase separation for large K is ex- 
pected since the ring term provides effective attraction be- 
tween bosons, see Appendix lB 3l and Ref.[l7'Sec. Vll. Also, in 
the large K limit, we can solve the ladder Hamiltonian exactly 
and find phase separation - this is shown in Appendix IB 21 
Nevertheless, for intermediate K the various phases discussed 
above are stabilized by the boson hopping. 



IV. CONCLUSIONS AND FUTURE DIRECTIONS 

An obvious virtue of searching for the fingerprint of pu- 
tative 2D gapless quantum phases in quasi- ID models is the 
relative numerical and analytic tractabihty of the latter An 
important requirement of this approach is that the 2D quantum 
phase has singular correlations along surfaces in momentum 
space. In this case the number of ID gapless modes present 
on an A^-leg ladder will necessarily scale linearly with N. 
While the particular fingerprint of the 2D quantum phase can 
already be present on the 2-leg ladder, as was the case for the 
DEL phase studied exhaustively in this paper, a compelling 
approach to 2D will certainly require analyzing models with 
larger TV. 

Already when iV = 3 exact diagonalization studies become 
problematic since finding the ground states of boson or spin 
models with more than 36 sites becomes virtually impossi- 
ble. Thus a 3-leg ladder of length L = 12 is already near 
the limit. DMRG is much more promising, and past studies 
on spin models reveal that convergence to the ground state is 
possible for lengths L of order of a 50 — 24 with iV = 3 — 6, 
respectively. 

Variational Monte Carlo for Gutzwiller-type wavefunctions 
constructed as products of determinants or Pfaffians as con- 
sidered in this paper is possible for significantly larger N, 
and equal time correlation functions can be computed with 
a fair degree of accuracy. But VMC is only as good as the 
variational wavefunction being employed, and the quality of 
the wavefunction is hard to assess in the absence of other 
more accurate methods. The ground-state energy of all phys- 
ically reasonable Hamiltonians is sensitive almost entirely to 
very short-range correlations, whereas for gapless quantum 
systems it is usually the longer range correlations which are 
necessary in order to distinguish between various competing 
phases. Quantum states in 2D with gapless critical surfaces 
in momentum space are somewhat of an exception since the 
location of such surfaces is presumably due to short-range cor- 



relations in the signs and amplitudes of the wavefunction. But 
the precise singular structure on such surfaces reflects longer 
range correlations in the wavefunction. 

Gauge theor}^ seems to offer one of the few analytic ap- 
proaches to access putative gapless spin-liquid type quantum 
phases in 2D that have no quasi-particle description. And for 
mean field states having a spinon Fermi sea, a controlled in- 
clusion of the gauge fluctuations is problematic, even in the 
simplest case of a compact U{1) gauge field. Accounting for 
the instanton (monopole) events which reflect the compact na- 
ture of the gauge field is exceedingly difficult. Even for a non- 
compact U{1) gauge field coupled to a Fermi sea, while there 
does exist an RPA-type approach, ''"''"''-^ it involves uncon- 
trolled and worrisome approximations. 

But for the ladder systems, gapless fermions (spinons) cou- 
pled to a ?7(1) gauge field is eminently tractable, although this 
fact has been exploited rather infrequently J2i22i2L22 In ID the 
gauge field has no transverse component and the bosonization 
method can be employed, leading to a gaussian effective field 
theory for the low energy excitations. This is another distinc- 
tive advantage of studying ladder systems relative to their 2D 
counterparts. 



A. DBL Phase on the N-Leg Ladder 

One extension of this paper involves studying the boson 
ring Hamiltonian on ladders with iV = 3, 4, 5 and larger. Al- 
ready for = 3 there are a number of new features that are 
not present on the 2-leg ladder First, one can consider either 
periodic or anti-periodic boundary conditions for the bosons 
in the rung-direction, and it is quite plausible that the phase 
diagram and the presence/absence of the DBL phase will de- 
pend on this choice. Second, while our search for a DBL 
phase with three gapless modes on the 2-leg ladder was unsuc- 
cessful due to instabilities of the DBL[2,2] variational state, it 
seems likely that for iV > 3 a DBL phase with more gapless 
modes than N should be accessible. This is interesting since 
in such a case at least two of the ID modes would have the 
same transverse momentum. 

For a general A^-leg ladder one can construct variational 
states of the form DBL[n, m] with integer n, m satisfying 
> ?7 > m > (with ji, m the number of partially filled 
bands for the di, d-z partons, respectively). The number of ID 
gapless modes in such a phase is ti + m — 1, so that DBL[3,2] 
for N = 3, say, would have four ID modes. If m = 1 the d2 
fermion has only one band partially occupied, and assuming 
the band has ky = 0, the corresponding variational wavefunc- 
tion vanishes whenever two or more bosons occupy a rung. 
In this case of no-double rung occupancy, the determinant for 
the c?2 fermion can be viewed as a (strictly) ID Jordan-Wigner 
string multiplying the di determinant. Whenever m > 1 this 
is no longer the case, and the d2 determinant would have a 
more subtle effect on the signs of the boson wavefunction. It 
would thus be desirable to access a DBL phase with m > 1. 
As detailed in Appendix I A 31 the gauge theory suggested an 
instability of the DBL[2,2] state for N ^ 2 due to a rather 
special "nesting" condition present for n — m ~ 2. For 
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n > m > 2 this will not generally be the case so we expect 
that such DBL phases should be more easily accessible for 
= 3, 4, .... It would also be interesting to use the DMRG 
to access some dynamical information about the DBL, which 
should be possible at least for = 2, 3 For N ^ oo the 
ladder model recovers the full symmetry of the 2D square lat- 
tice, and in future work it would be desirable to study how to 
approach this limit. 



B. SU(2) Invariant Spin Models with Ring Exchange 

It would be most interesting to search for possible 2D 
spin liquid phases with singular surfaces in momentum space 
in models possessing SU(2) spin symmetry. Recently, sev- 
eral author o'^i'^i'^ have suggested that the putative spin liq- 
uid observed in k-(ET)2Cu2(CN)3 is of such character, and 
have proposed that a mean field state with a Fermi surface of 
spinons is an appropriate starting point. Variational energet- 
ics on the corresponding Gutzwiller projected spinon Fermi 
sea have been performed on the triangular lattice Heisenberg 
antiferromagnet with a four-site cyclic ring exchange term, a 
model argued to be relevant for a Mott insulator such as k- 
(ET)2Cu2(CN)3 that has a comparatively small charge gap. 

It should be possible to study the Heisenberg plus ring ex- 
change spin Hamiltonian on the triangular strip to search for 
signatures of the proposed parent 2D spin liquid. In the ab- 
sence of ring exchange this model is equivalent to the ID 
Heisenberg model with both first and second neighbor inter- 
actions, Ji, J2, and has been explored using DMRG in ear- 
lier work.i^^i^ The phase diagram as a function of J2 / Ji ap- 
pears to have two phases: a Bethe-chain phase stable for small 
J2 / Ji with one gapless mode and a fully gapped dimerized 
state for intermediate and large J2/ Ji ; the gap decreases ex- 
ponentially for large J2 / Ji where we have nearly decoupled 
legs ( J2) of the triangular ladder (with interchain Ji). 

With a 4-site ring exchange term (coupling K) the triangu- 
lar strip model is characterized by two dimensionless param- 
eters, J2/J1, A'/Ji. Are there any additional phases in this 
phase diagram? The spin liquid phase on the triangular strip 
that descends from the (putative) 2D spinon Fermi surface 
state should have 3 gapless modes (2 for spin, 2 for the two 
legs, minus one from the gauge constraint) and corresponding 
incommensurate spin correlations. Variational Monte Carlo 
together with a gauge theory analysis can provide a detailed 
characterization of this state, which could then be compared 
with DMRG, precisely as was done in this paper for the boson 
ring model. If present, one expects the spin liquid phase to 
appear for intermediate values of both J2/J1 and A / Ji. 

In addition, one might study a half-filled Hubbard model 
on the triangular strip, which should exhibit a metal insulator 
transition at intermediate coupling, U/t. Of particular interest 
is the nature of the quantum state just on the insulating side of 
the Mott transition. If the Mott transition is weakly first order, 
there will be substantial charge fluctuations and ring exchange 
interactions in this part of the phase diagram, and possibly a 
spin liquid phase, in addition to phases that have already been 
identified (see Refs. [3lll32ll33l and references therein). 



C. Itinerant Electrons 

What are the prospects of using ladders to approach non- 
Fermi liquid phases of 2D itinerant electrons? One complica- 
tion is as follows. Imagine a weak coupling Hubbard model 
on the square lattice (at densities well away from special com- 
mensurate values) which has a conventional 2D Fermi liq- 
uid ground state. On an A^-leg ladder the free Fermi points 
present for [/ = would be converted into Luttinger liq- 
uids. The jump discontinuity in the momentum distribution 
function at each Fermi point would be lost, but singularities 
would remain, characterized by some Luttinger liquid expo- 
nents. Moreover, the location of the singularities would still 
satisfy the Luttinger sum rule. Now imagine some strong cou- 
pling 2D electron Hamiltonian with a non-Fermi liquid phase 
that has a residual Luttinger Fermi surface but with Z = 0, 
analogous to a ID Luttinger liquid. The corresponding A^- 
leg ladder descendant would be qualitatively indistinguish- 
able from the phase of the weak coupling Hamiltonian. Con- 
versely, the presence of a Luttinger-satisfying Fermi surface 
on an iV-leg ladder would not enable one to distinguish be- 
tween the two different 2D phases, one a Fermi liquid and the 
other not. 

On the other hand, the A^-leg ladder descendant of a 2D 
non-Fermi liquid ground state with momentum space singu- 
larities that violate the Luttinger sum rule (surfaces with the 
"wrong" volume, or perhaps even arcs) would have quali- 
tatively distinct signatures. Several recent paper a'^-'^ have 
proposed such 2D non-Fermi liquid phases, and it would be 
extremely interesting to find evidence for their ladder de- 
scendants. The wavefunction for one of these phases was 
constructed^^ by taking the product of a free fermion deter- 
minant and the DBL wavefunction in Eq. This phase, 
which inherits the d-wave sign structure from the bosonic 
DBL wavefunction, was called a d-wave metal and would ex- 
hibit distinctive signatures if present on a ladder system. A 
possible Hamiltonian that might possess the d-wave metal can 
be expressed by adding an "itinerant-electron-ring" term to 
the usual square lattice t — J model: H = Ht-j + H^l^^g 
with, 

^^rtg = ^el ^(5^,352,4 + H.C.), (21) 
□ 

where the sites 1, 2, 3, 4 are taken to run clockwise around an 
elementary square plaquettes and the summation runs over all 
packets. Here we have defined an electron singlet creation 
operator: 

S^,,^{cl^cl-clc]^)/V2. (22) 

The ring term rotates a singlet on the 2 — 4 diagonal into a 
singlet on the 1 — 3 diagonal. It would be interesting to study 
this or other such strong coupling Hamiltonians on the iV-leg 
ladder. 
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D. Conclusions 

To conclude, in this paper we have initiated a study of 
A^-leg ladder systems that are descendants of candidate 2D 
quantum phases with low-energy excitations residing on sin- 
gular surfaces. Motivated by one such proposal for the so- 
called DBL phases of un-condensed itinerant bosons,'^ we 
have studied the 2-leg model with frustrating ring exchanges 
using exact numerical approaches. We have indeed found the 
DBL[2,1] ladder version arising prominently in the phase di- 
agram. We have also searched for the DBL[2,2] version but 
from the gauge theory analysis concluded that it is likely un- 
stable due to special kinematics on the 2-leg ladder and in- 
stead gives rise to a boson-paired phase with extended s-wave 
character in the pairs. This paired phase is in fact realized 
in our model, while the DBL theory gives us tools to under- 
stand its properties. While our focus has been on the DBL 
ideas, we have explored the full phase diagram of the specific 
ring model in fair detail and found also other strong-coupling 
phases such as the above s-wave paired phase and the d-wave 
paired phase. The latter is not directly accessible from the 
DBL theory but is characteristic of the binding tendencies in 
our ring terms, which eventually cause the bosons to phase- 
separate for large ring exchanges. We hope to pursue similar 
ideas in the > 3 ladders and also in SU(2) spin and itinerant 
electronic models with ring exchanges, which are particularly 
exciting. 
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APPENDIX A: GAUGE THEORY DESCRIPTION AND 
SOLUTION BY BOSONIZATION 

A faithful formulation of the physical system in the slave 
particle approach Eq. ([T]i is a compact U(l) lattice gauge the- 
ory. We set this up on the 2-leg ladder as follows. Denote the 
vector potential components on the links of the two chains as 
ai,x and au^x, and on the rungs as Uy. In the Euclidean path 
integral, denote the temporal components associated with the 
sites on the two chains as a/ and a// The action for the 
gauge field contains — cos{aj,x^aii.x+'^xay), — cos(a7^^ — 
aii.T + ^ray), and - X]/=/,// cos(Vra;^j; - Vxai^r)- For 
each / = /, //, the last term gives a 1h-1D gauge action, and if 
the cosine is interpreted as a Villain cosine, we can replace it 
with (Vrft;,!; — Va;a;^7-)^ and treat a; as non-compact with no 
change in the results. Choosing the gauge ay = 0, the combi- 
nation a J — aji is massive, while the combination aj + ajj is 
massless. 



Having described the gauge sector, we now consider the 
fermions di and d2 with the mean field bands as in Fig. |2] 
and take the continuum limit using the bonding da^ and 
antibonding d'^'' fields near the corresponding Fermi points 
P — R/L = +/—■ These fields couple to a = (a/ + a//)/2 
as to the usual gauge field, while the massive combination 
0.1 — ail can be integrated out. The continuum Hamiltonian 
density is 



^kinetic ~t~ 



^kinetic 



a,ky,P 



Pv 



Fa "'aP 



(Al) 

{-tdx - e^ax)d'^!:'p\{A2) 



where ei = 1 and e2 = — 1 are the gauge charges of the di and 
d2 respectively. In the DBL[2,2] case, there are both bonding 
and anti-bonding fields present for each species di and d2- 



The Fermi momenta satisfy k 



(0) 



27rp, where p is the 



original boson density per site. In the DBL[2,1] case, the ^2 
fermions have only bonding fields. 

The allowed interactions contain general density-density 
terms, 

/i*-/,L_ = ^ [Bb^b'PbPPb'P + Cb^'PbPPb'-p] , (A3) 



if erm 



b,b'.P 



where 6, b' sum over all bands of all species, and also the fol- 
lowing terms. 



,^{11) - 

"-if erm ~ 



^(0)t ,{0)t .W ,W , TT 



F 



G 



j:df»d 



(")t^(0)t 

2-P"l-P"2P 



(o)tjWt AO) 



-pdl 



P"2P 



H.c. 



H.c. 



(A4) 
(A5) 
(A6) 



In the DBL[2,1] case, all terms that contain ^2 are absent. 

The above QED2-like theory can be analyzed perturba- 
tively in the matter-gauge coupling, e.g., in a systematic 1 /N 
expansion, as was done in Ref. [191 for the spinon-gauge treat- 
ment of the ID Heisenberg spin chain. From such studies, 
we often use the following "Amperean" rule of thumb: Gauge 
interactions modify in a singular way processes that involve 
fermion fields with oppositely oriented group velocities. If the 
corresponding gauge currents are parallel (antiparallel), such 
processes are enhanced (suppressed), which originates from 
the attraction (repulsion) of currents in electromagnetism. As 
an example, d^^d^^L oppositely charged particles moving 
in opposite directions and producing parallel gauge currents, 
so this bilinear is expected to be enhanced compared to the 
mean field. We will see this explicitly in a solution of the ID 
gauge theory by bosonization,'^'-^ which we pursue instead 
of the perturbative treatment. We also caution here and will 
see below that the different Fermi velocities and the general 
allowed density-density interactions complicate the analysis 
significantly and can independently affect power laws, in ad- 
dition to the above Amperean rule for the gauge interaction 
effects. 
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To bosonizer^i^^^ we write 



(A7) 



Here 



b^a"^ and e^t" 



are the conjugate phase and phonon 



fields for each band, while t]^ are Klein factors, which 
we take to be Hermitian operators that commute with the 
bosonic fields and anticommute among themselves. The ki- 
netic Hamiltonian density becomes 



^kinetic 



\ " "Fa 



cx,ky 



2tt 



(A8) 



The density-density interactions Eq. iA3i lead to generic 
terms of the form {dx4'b — ebax){dx4>b' — e.b'ax) and 
dxdbdx(^b', which are strictly marginal and affect power laws. 
The other interactions Eq. ( IA4l i produce cosines and are writ- 
ten separately for the DBL[2,1] and DBL[2,2] below. 

We proceed in the Euclidean path integral, choose the gauge 
Ux = and integrate out the ar- This renders the field 
6c ^ J2a k ^q'^q""'' massive and essentially pins it to zero. 
This is the suppression of the charge fluctuations in the gauge 
theory and realizes the microscopic constraints d\di = d\d2 
at long wavelengths. We then have two gapless modes left in 
the DBL[2,1] and three modes in the DBL[2,2], if we can as- 



are irrelevant 



sume that the cosines corresponding to nifcrm 
(see below). 

To characterize the phases of the original hard core bosons, 
we examine the single boson correlations, the boson density 
and current correlations, and the pair boson correlations. 

The microscopic boson field is written as 



b{x,y) = di{x,y)d2{x,y) 



E 



(A9) 



kly,k2y,Pl,P2 

j{Pt k'^^^y ' +_P2^>|'' ' )xrik,y+k2y )vAkly)Ak2y ) 



which we then express in terms of the continuum bosonized 
fields (see separate treatments of the DBL[2,1] and DBL[2,2] 
below). For oppositely moving Pi and P2, the 9 part has a 
non-zero projection onto 9c, and gapping out the 9c by the 
gauge fluctuations enhances the corresponding contribution, 
in agreement with the Amperean rule applied to the oppositely 
charged di and d2- 

The fermion densities are 



Pd^{x,y) = dl{x,y)dc{x,y) 



(MO) 



P,P',ky 

p,P' 



the rungs (the definition below uses the gauge choice ay 
made early on): 



J±,a(a;) - - [4(a:,l)dja;,2) -4(x,2K(x,l)](All) 



P.P' 



i(Pfc^ 



The continuum expansion resembles that for the density at 
qy ~ IT, but with opposite signs between df^d^^p, and 

'^q^p''^q'-P' '■^i^ ^ill h^lP distinguish different phases in 
the analysis of instabilities, Secs. lA 2llA3l Bosonized expres- 
sions for the contributing bilinears will be given separately in 
the DBL[2,1] and DBL[2,2] cases. Here we only note that 
when the constituent particle and hole are on the opposite 
sides, we again have non-zero projection onto 9c, so the corre- 
sponding contribution is enhanced by the gauge fluctuations, 
in agreement with the Amperean rule. 

Finally, we consider the boson pair operator 
b{xi,yi)b{x2,y2) and expand in terms of the continuum 
fermion fields. Assuming xi and X2 are nearby, we character- 
ize the pair by its center of mass coordinate X — {xi +X2)/'2 
along the chains (this becomes the argument of the continuum 
fields) and also by the internal pair structure. To this end, we 
collect different microscopic contributions that give rise to 
each term dia{X)dib{X)d2c{X)d2d{X), where a, b, c, d are 
combined band and Fermi point indices. Using short-hands 
Diair) = diaexp[ikia ■ r], etc., then b{xi,yi)b{x2,y2) 
contains 

- [Diaixi,yi)Dib{x2,y2) - Dia{x2,y2)Dib{xi,yi)] 

X [D2c{xi,yi)D2d{x2,y2) - D2cix2,y2)D2d{xi,yi)] 
= dia{X)duiX)d2c{X)d2d{X) 

where £^ ~ xi — X2 is the relative coordinate along the chains. 
The "pair wavefunction" is 

$la,lfc,2c,2rf(yi,y2,e) =e'''=^-+'^"''+'^^-+''^='^''^'' (A12) 
. /kia-kifc \ . /k2c-k2d \ 

X 4 sm I • p sm ' P j 



with Y = 2£i±y2. and p = j/i — 2/2)- In the above, we have 
chosen to characterize the pair by the momentum along the x- 
axis, but not by the momentum along the y-axis. On the 2-leg 
ladder, it is easier to visualize pairs by keeping track of both 
yi and j/2- 

A few more words about this characterization on the 2-leg 
ladder. Consider first a 2D setting, where we would write a 
contribution to a boson-pair 6(ri)6(r2) by say 



(A13) 



The two lines in the expansion separate qy = and qy = tt 
parts, and all contributions with the same non-zero momentum 
are grouped together (assuming generic distinct bonding and 
antibonding bands). We also consider the particle current on 



where O is some slowly varying operator of the center of 
mass coordinate R, while $(p) is the internal pair function 
in the relative coordinate p. Specializing to a square lattice, 
we could then distinguish s-wave and d-wave (more precisely. 
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dxy) pairing by looking at the signs of $(p = i ± y). On the 
2-leg ladder, however, separating out the transverse momen- 
tum Qy as in 2D mixes things a httle. Specifically, a rotation 
of an a; + 2/ pair into an x — y pair can be also achieved by 
a "translation" y y + 1 mod 2. To avoid this, we instead 
define the pair-function by writing 



021cg(X)e*«-^$21cg(^,yi,2/2), 



(A14) 



where X and <^ are the x-components of R and p respectively. 
We then characterize the pair wavefunction $2iog by the sym- 
metry under interchanging the two legs of the ladder {yi.2 
?/i^2 + 1 mod 2). When $2icg is even, we call it s-wave, since 
then the amplitudes for a pair sitting on a +45° diagonal and 
a pair on a —45° diagonal are the same. When <f>2iog is odd, 
we call it d-wave, since then the amplitudes for the two diago- 
nals are opposite. The names "s-wave" and "d-wave" are used 
in anticipation that it is such local pictures of the pairs on the 
plaquettes that survive when we build towards 2D. Returning 
to the DBL construction Eq. ( IA12I ). we get s-wave or d-wave 
character depending whether kiay + fcitj, + k2cy + k2dy is even 
or odd multiple of tt. 

We now specialize to the DBL[2,1] and DBL[2,2] cases in 
turn. 



involved biUnears by their momentum: 

(0,0): 



/ . "ip "IP ^ 



(24"i'\ 0) : d^^fd^^l^ = e-'^'^+'^^'- 
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(0,0) : 
(2fcf2,0) 



Here we include the Klein factors, but do not show combi- 
nations that can be obtained from the above by reversing all 
momenta. 

When we include gauge fluctuations, 6c becomes massive, 

locking together df^ + and 02- Integrating out ((/)c, ^'c), 
we obtain a generic harmonic liquid theory in terms of two 
coupled modes ((/)p, Qp) and (0_ , BJ). Of the cosine interac- 
tions Eq. (IA4l i. we have only 



M/erm = 2£^1 COs(40_) 



(A15) 



1. DBL[2,1] 



a(0) 



Upon bosonization, we start with three modes y^-^ 

((\)\^ , ^'i'^'' ), and ((/)2 , 6*2) (we drop the bonding band label for 
the d2 fermions). To proceed formally, we change to new 
canonical variables via 
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202), 



(0) 



^>l^'-'^2), 



1 



2//)(0) 



202 



02); 



9W. 



The microscopic boson field 6(r) has the following contri- 
butions listed by their momentum (cf. Eq. 



{k^p^^ +kF2,ky) : df^'>d2B. 



^F2 



ky) 



Iip "2P 



□ ±0-)gi(' 



In each equation, the top and bottom signs are for ky = and 
TT respectively. We omit the Klein factors for simplicity since 
we will not need them here. Also, we do not show combi- 
nations that can be obtained from the above by reversing all 
momenta. 

The microscopic boson density has contributions from both 
the di and ^2 particles via Eq. jAlOl i, while the rung current 
has contributions only from the di via Eq. ( lAllI ): we list all 



Assuming this is irrelevant, we have a stable phase with two 
gapless modes. The potential instability due to this term is 
considered in Sec. lA 2l below. 

We do not write explicitly the full DBL[2,1] theory in terms 
of the {(j)p, Op) and ((/)_ , 6'_ ) since even without the additional 
interactions from Eq. ( |A3t we have a general coupled har- 
monic system. We can still make some observations about 
the scaling dimensions of the contributions to the boson and 
the boson density or current operators. In the initial fermionic 
mean field, the bands d^i^'"\ d2 are decoupled and the scal- 
ing dimensions of all contributions to the operators h and pf, 
are equal to 1 . The gauge fluctuations effectively set 9c ~ Q 
everywhere. 

Consider first the contributions to h. For arbitrary couplings 
in the full harmonic theory, we can argue that the scaling di- 
mensions of the operators d\p d2B. are at least 1, i.e. larger 
than the mean field value (these have no explicit Amperean 
enhancement). On the other hand, setting 9c = 0, the scaling 

(k ) 

dimensions of the operators d\p d2L can be lower than the 
mean field in accord with the Amperean rule; without further 
information, we cannot say much - in general, these scaling 
dimensions can be as low as zero. 

Consider now the contributions to the boson density pf,. 
The zero momentum contribution has scaling dimension 1. 
We can argue for arbitrary couplings in the theory that the 
{k'^p\ — k'^pl , tt) contribution has scaling dimension larger than 
the mean field value 1. On the other hand, the (2A:^j"\ 0), 



{k^p\ + k'^p-l, tt), and {2kp2, 0) contributions can have smaller 



scaling dimensions upon setting 



0, in accord with the 
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Amperean enhancement rule. We do not have general bounds 
on the scaling dimensions in the coupled two-mode system, 
although clearly these are not all independent. Furthermore, 
in all of the preceding discussion, even if the gauge interaction 
acts to enhance some correlation, the short-range interactions 
in ID can act to suppress the correlation, so the ultimate fate 
is not clear 

As an illustration, we list all scaling dimensions in the case 
when the {<f)p,6p) and {(p_,9_) modes are decoupled and 
characterized by the Luttinger parameters gp and g_ respec- 
tively [our convention is that g enters as g{dx4>Y + 
in the action]; 



chains: 

: 

P{kpl — kp^) 



d^iR d'll' d2Rd2L 



'■IL 



,i(20p±20_) . 



"lP"lP"2fl«2L 

d^ipd^i}pd2Rd2L 



As explained earlier, we characterize each contribution by the 
internal pair wavefunction Eq. ( IA12I ). Again we have omitted 
the Klein factors but can easily restore them when needed. 
For the zero momentum contribution we have 



^lR{ky),lL{ky),2R.2L{yi,y2,0 — 

Ae'^y^y^+y^-^ sin(4^''^0 sin(fcF20 



(A16) 



^K,(24'=-\0)] 
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1 1 
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9p 
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This example is loosely motivated by the observation in 
Fig.[8]that the boson n(q) has roughly similar singularities at 

{k^p\ — fcF2, 0) and {k^pl — fci?2, ti")- Also, increasing the ring 
term K and decreasing the interchain hopping J± drives the 
mean field bonding and antibonding bands to be more similar, 
suggesting such approximate decoupling of the two modes. 
Continuing with this illustration, the boson correlators at 

fcF2, ky) are most enhanced when gp ^ 2 and g^ ^ 2 
- 1 /2; with such couplings we also 



giving A[6^^(.,,_^^^^^^^ 
find mean-field-like A[p ] 



1, while A[p, 



'b,{2kf 



,0)] 



2. These 



numbers are qualitatively similar to the singularities de- 
tectable by eye in the density structure factor in Fig. |8] where 
we see some signatures at (2/s^/^-', 0), (fc^^^ ± k^p^,TT), 
and no visible signature at (2fci?2,0). A more careful look 
at the DMRG data shows that the density singularities at 
{2kp(^\o) are slightly stronger than the mean field, while 

the ones at {kpl ± k^pl^-n) are weaker; this would suggest 
g- somewhat smaller than 2, which would also be consis- 
tent with the assumed irrelevance of the cosine interaction 
Eq. (IA15b . However the DMRG boson correlator singular- 
ity has estimated scaling dimension somewhat smaller than 
the maximum 1/2 in this illustration, so one should not take 
the above too literally. Nevertheless, this example gives some 
sense to the strengths of the observed singularities in Fig.|8] 

Finally, consider a boson pair operator; we list contribut- 
ing four-fermion combinations by their momentum along the 



This remains unchanged if we simultaneously change the y 
coordinates of both bosons in the pair, e.g., if we move a pair 
that lies entirely in chain 1 to chain 2 or if we turn a +45° 
diagonal pair into a —45° diagonal pair; we call this s-wave. 

2-Kp we 



For the contribution at momentum k 



have 
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(A17) 

- y2) 



uio) _ uy^i 



This changes sign if we simultaneously change both y coordi- 
nates in the pair When kpl and k^p^ approach each other, the 
pairs straddling the two chains have larger amplitude. Since 
the amplitude changes sign when we turn a +45° diagonal 
pair into a —45° diagonal pair, the internal structure of the pair 
has some d-wave character The contribution at kpl — k^p^ can 
be similarly characterized, but is omitted here. 

In the free fermion mean field all contributions have scal- 
ing dimension 2. Beyond the mean field, the d-wave contribu- 



tion at kpl + fc^/ is potentially enhanced by setting 9c = 0; 
still, its scaling dimension is at least 1, while we cannot say 
much about the scaling dimension of the s-wave contribution 
at zero momentum. In the above illustration with decoupled 
<j>p and 0_ modes, the d-wave pair has scaling dimension 
1/gp + .gp/4 = 1 if we set gp = 2, while the s-wave pair has 
scaling dimension is 1/gp + l/.g- = 1 if we also set g^ = 2. 
The latter would be larger if we take smaller g-, in agreement 
with the observed dominance of the d-wave-like pair correla- 
tions in the DMRG Fig.[TO]oscillating at wavevector 27rp. 



2. Possible Instability of the DBL[2,1] 

Let us ask what phase is obtained starting from the 
DBL[2,1] if the interaction Eq. ( IA15I ) is relevant [the Am- 
perean rule does not point one way or another, but additional 
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interactions from Eq. ( IA3I ) can make this happen]. Then </)_ 
is pinned, with two distinct cases depending on the sign of Ei 
considered below, while 6'_ fluctuates wildly. All contribu- 
tions to b contain 9^ in the exponent, so the single boson cor- 
relation function decays exponentially - there is a charge gap 
to single boson excitations. The {(j>p,9p) mode still remains 
gapless, and we have a single-mode harmonic fluid described 
by the parameter gp. 

The boson density shows power law correlations at zero 
momentum with scaling dimension 1 . The remaining particle- 
hole bilinears with power law correlations are at wavevector 
{kpl + k^pl, tt) = (27r/9, tt) (scaling dimension gp/4) and at 
(2fci?2,0) = (47rp, 0) (scaling dimension gp). The q = 
and q = {2kp2,Q) parts do not distinguish qualitatively be- 
tween the two possibilities Ei > {) and Ei < 0, while the 



{k 
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k^Fi ' P^'"'^ Aoe^^ distinguish. Explicitly, the di 



particle density is 

pa^{x,y) = -2irif^ri'^;'^ siii(2(/._)e'^^e'('=^"^+'^^"')"+*"^ + H.c., 

where we have omitted the zero momentum component and 
parts with exponentially vanishing correlations. On the other 
hand, the di particle rung current is 



Consider now two cases: 
El > 0: It follows that 24>^ 



J mod TT. In this case p^i 



shows power law correlations at (kpl + k'^p^, tt), while j^^i 
correlations are absent. 

El < 0; It follows that 2(/)_ = mod tt. In this case p^i 
correlations at the above wavevector are absent, while j±j 

shows power law at wavevector kpl + k^p^ along the chains. 

Let us finally consider the pair-boson correlations, which 
are also power law. Putting together all contributions to 

b{xi,yi)b{x2,y2),v^e have 

e'^"*''[e'2'^-$lK(0),lL(0),2i?^,2L + e~'2"^-$ii?^(^),iLW,2i?^,2L] 



where the appropriate ^{yi ,y2,£,) are given in Eqs. ( IA16I ) and 
( IA17b . The pairs at zero momentum have dominant correla- 
tions with scaling dimension 1 /g^; these have s-wave char- 
acter as far as rotating diagonal bonds is concerned, but the 
amplitude details depend on the sign of Ei. The pairs at 
kpl+k'^p^ have subdominant correlations with scaling dimen- 
sion l/gp + (7p/4; these have d-wave character for rotating the 
diagonal bonds and the details do not depend much on the sign 
of£:i. 

Summarizing, the resulting phase has a gap to single boson 
excitations, but has power law correlations for particle-hole 
and particle-particle composites. The dominant particle-hole 
composite has scaling dimension (?p/4 and either contributes 
to the density at wavevector (27rp, tt) if Ei > or to the rung 
current at 27rp if Ei < 0. The dominant particle-particle com- 
posite has scaling dimension l/gp and represents s-wave pair- 
ing at zero momentum. Loosely speaking, we can describe 



this phase by saying that the bosons form s-wave-like pairs 
and these molecules in turn form a Luttinger liquid. The cur- 
rent and/or density fluctuations occur at a wavelength equal to 
the mean inter-pair spacing along the ladder, as expected for 
a ID Luttinger liquid (of pairs). The form of the particle-hole 
fluctuations across the rungs presumably reflects the internal 
pair structure in the two cases. 



3. DBL[2,2] and Instability Towards S-Wave Pairing 



In the DBL[2,2] case, we start with four bands, 
d^i^^\d^^^\ and have four modes upon bosonization 
Eq. ( IA7I ). It is convenient to introduce the following canonical 
variables: 



1 

2' 
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2^ 
1 

^(ri 
1 
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^'i'+ 






,(0)„ 




4"' 


+ 4" 



(A18) 

(A19) 
(A20) 
(A21) 



with the same transformation for the 9 variables. The cosine 
interactions Eq. ( IA4b are 

4Srm ^ 2£;i cos(2</)+ + 20_) + 2E2 cos(20+ - 20_) 
-4fFcos(2(/i+)cos(26'c) +4fG'cos(2(/)_)cos(26'c) , 

where f = Vi^^ Vi^"^ V2^'' ''12''' ■ 

Proceeding with the analysis as in the DBL[2,1] case, upon 
integrating out the gauge field, 9c becomes massive. If by 



tweaking the strictly marginal interactions h 



ill) 

Aferm 



(I) 

Aferm 



we could 



terms irrelevant, we would end up with a 



render the h 

phase with three gapless modes. We call this possible phase 
DBL[2,2], and it can be analyzed similarly to DBL[2,1]. For 
example, the mean field boson correlation in this phase would 
read 



cos[(fc(^r'-4'r')x] 



- cos[(A: 



Fl 



The first term in the brackets comes from <i[p«''4-^|,' and 
is expected to be enhanced by Amperean attraction (in the 
bosonization, this term contains 9c in the exponent and is 
potentially enhanced upon setting 9c = 0). On the other 
hand, the second term in the brackets has no Amperean en- 
hancement. As emphasized earlier, besides the gauge fluctu- 
ations crudely captured by the Amperean rule, other interac- 
tions can also change the scaling dimensions. Interestingly, 
and this has important consequences explored below, d^^^d'^l 
and d^^^d^il carry the same momentum {k]pl — kpl,0) = 



(^F2 ~ ^f7'0), and also d'^^'^d'^j^ and diX^d2fl carry the 



(0) j(^) 



MAO) 
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same momentum (/c^] — k'^p^jTr) ~ {k 
fact, Josephson-like coupling between the first two "bosonic 
modes" is allowed and gives rise to the F interaction in 
Eq. (IA4l i. while Josephson coupling between the last two 
modes gives rise to the G interaction. If the gauge physics 
is the dominant determining factor in the DBL[2,2] phase, 
we expect the state to be strongly unstable to these inter- 
actions, and we pursue this scenario below. On the other 
hand, if the DBL[2,2] were stabilized against the F and G by 
some density-density interactions dominating over the gauge 
physics, which is possible and interesting, we would have lit- 
tle intuition from the gauge theory perspective, and we do not 
pursue this possibility further. 

If the gauge physics dominates, it appears very likely that 
the F and G terms are relevant, as can be seen from their 
bosonized expressions upon setting 9c ~ (in accord with 
the Amperean rule). Let us explore the resulting phase when 
F and G flow to large values and pin the fields and . 
Four different possibilities depending on the signs of F and G 
are discussed below. In such a phase, the conjugate fields 
and 0_ fluctuate strongly and any operator containing these 
in the exponent will have only short-range correlations. Thus, 
we conclude that the boson correlations decay exponentially, 
i.e., there is a gap to single boson excitations. However, there 
is still one gapless mode {(j)p, Op), which we characterize with 
the Luttinger parameter gp. 

Consider fermion bilinears composed of a particle and 
a hole of the same species that enter density and current, 
Eqs. (lAlOIAlll i. Other than the <? = densities d^^^^'^ d'^!""^ 
the following bilinears survive when 0± get pinned: 



JO) 
F2 



In 



,(O)„(7r)„-i(0++e„0_)„-,:P(ep+e„e,) 



^p^'^a-p =C^reH*'++^°*'-'e-'^^''''+^°'^=^ (A23) 



M„{0) J{4'++e^<P-) „-iP{e,+e^e^} 



The particle density Eq. (lAlOb is, upon dropping the zero mo- 
mentum part and setting 6c = 0, 

, (x, y) ^ -2zr;i")r,i^) sin(,^+ + e„0„) (A24) 



The particle current on the rungs Eq. (lAllb is 



(A25) 



We now discuss different phases that can arise depending 
on the signs of the couplings F and G. We set F = 1, from 

which it follows (?]2"''^2^'') = ~{v'i'\i^'') (the physical re- 
sults are independent of the choice of F). There are four cases: 
F > 0,G > : It foUows that 0+ = mod tt, 0_ = f 
mod TT. In this case pd^ = Pd2 show power law correlations 
at (27rp, tt) with the scaling dimension .gp/4, while j±_a cor- 
relations are absent. This is a natural phase coming from the 
microscopic gauge theory since the fluctuations of the di and 
d2 densities are in sync; also, when we integrate out the mas- 
sive gauge field aj — ajj early in the derivation of the contin- 
uum theory, we in fact generate such F = G > couplings. 



Foretelling the analysis of the pairing correlations below, we 
propose that this is the s-wave-paired phase observed in the 
DMRG, Sec. lmCl 

F < O.G < 0: It foflows that = ^ mod tt, 0_ = 
mod TT. In this case = —pd2 show power law corre- 
lations at {2ttp,tt), while j± ci correlations are absent. The 
fluctuations of the di and d-z densities are out of sync. This 
phase with spatial modulation of the gauge charge is less nat- 
ural coming from the lattice theory. It is in principle allowed 
if the bare couplings are finite, since then such modulation 
costs only finite energy density that can be offset by some 
other short-range interactions; but the cost is large if the mi- 
croscopic theory is at strong coupling as is usually the case 
in the slave particle treatments. For example, only the phase 
with pd-^ ~ Pd2 can be realized by our wavefunctions. 

F < 0, G > 0: It foflows that 0+ = 0_ = | mod tt. In 
this case pd^ correlations are absent, while = j±,2 show 
power low correlations at wavevector 2ttp along the chains 
with scaling dimension .gp/4. The di and d2 particle currents 
are in sync, so this phase is natural in the gauge theory; for 
the original bosons, it would have enhanced current-current 
correlations but not density correlations. However, we have 
not observed such a phase in the ring model. 

F > 0, G < 0: It follows that 0+ mod tt. In 

this case pd^ correlations are absent, while jj^.i = —j±.2 
show power low correlations at wavevector 2np along the 
chains. Since the currents are out of sync, this phase is less 
natural coming from the microscopic gauge theory. 

Let us consider boson pairs. First, for each a = 1,2 
there are four dad^ continuum field combinations that survive 
when the are pinned: 



"'aP"-aP — 



^iO)^iJ)^^i<P,+e.4>.)^^Pi0,+e.e^)^ (A26) 



pi[0p+e„0c±(0++eo<^-)] 



(A27) 



where in the last line the top/bottom sign corresponds to 
ky = O/tt. Constructing boson pairs via didid2d2, the field 
(pc disappears, and we also set 6c = 0. Of the sixteen com- 
binations, six carry zero momentum and aU have the largest 
scaling dimension 1 / g^: 



^(0) .(0) M 

"iP^lP "2-P"2-P 

ji^y) ji^y) ji^y) ji^y) z{2(^p±20+) 

"ifl "iL "'2P "2L — ^ ' 

jiky) i{ky) ,iky+TT) ,{ky+Tl) _ i(20p±20_ 
"iP "iL "'2P "2L — 



(A28) 

(A29) 
(A30) 



The second and third lines will likely have smaller amplitudes, 
since cp^^ and (/)_ still fluctuate a little about the pinned values. 

We now examine the internal pairing structures. Corre- 
sponding to the first line above, we have: 

*J'lP(0),l^"(7r),2-P(0),2-P(7i-)(yi, 2/2,0 =(A31) 

= -ASy.^y^ sin ( ^£L_^^ 1 sin ( ) 



-45, 



'yi¥=v: 



cos 



'l(o) _ -(^) \ /i.(o) _ . w ^ 

"■PI '^Pl c j (jQg ^P2 "'P2 c 
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The Sy-^^y^ piece describes a pair formed by bosons in the 
same chain, while 5y^^y^ is a pair straddling the two chains. 

When say fc^^^ « which becomes more accurate as 

we decrease J^, the amplitudes for the latter are signifi- 
cantly stronger, i.e., pairs are predominantly straddling the 
two chains. 

Corresponding to the other zero-momentum combinations, 
we have 

*lH(fe„).lL(fcJ,2_R(fe;),2L(fe;)(2/l,y2,6 = (A32) 



4e 



i(ky+k' )(yi+y2) , 



If ky = k'y, this is independent of whether the two sites are on 
the same or different chains. On the other hand, if ky 7^ ky, 
this has opposite signs if the two sites are on the same versus 
different chains. 

All of the above pair-functions are even if both particles are 
moved perpendicular to the chains. In particular, +45° and 
—45° diagonal pairs have the same amplitudes. This is what 
we call s-wave pairing. Putting together all zero momentum 
contributions to b{xi,yi)b{x2,y2), we have 



2r$iij,(o),li?,(7r),2L(0),2L(7r) 

e'^'^+$li?,(0),li(0),2i?.(0),2L(0) 
e 1*li?.(7r),li(7r),2i?.(7r),2L(7r) 



-i20_ 



lfl(0),lL(0),2i?.(7r),2L(7r) 
*lfl(7r),lL(7r),2fl(0),2L(0) 



(A33) 
(A34) 
(A35) 
(A36) 

(A37) 



where F = rj^^ r][^^ t]^'' 1)2^^ and eventually drops out since 
both r cos(2(/)+) and rcos(2(/)_) are determined uniquely by 
the signs of the F and G couplings: re**2'^+ = sign(F), 
Pg±i20_ _ _sign(G'). Note also that the contributions in the 
last four lines will have smaller amplitudes since (/)-|_ and 0_ 
fluctuate a little around the pinned values. So the numerically 
largest contribution is given by the first line, while the smaller 
contributions from the last four lines will add to it with relative 
signs that depend on the signs of the F and G. 

Summarizing, the DBL[2,2] is unstable towards a boson- 
paired phase with s-wave pairs carrying zero momentum. The 
resulting phase is roughly similar to the s-wave paired phase 
discussed as a possible instability of the DBL[2,1] in the pre- 
ceding section, but some details are different. For example, 
coming out of the DBL[2,1] the amplitude for rung pairs van- 
ishes, while out of the DBL[2,2] the rung pairs have large am- 
plitudes comparable to those of the diagonal pairs. The latter 
is more similar to what the DMRG finds in the s-wave paired 
phase in the ring model. Sec. IIIICI The prediction of dom- 
inant density correlations at wavevector {2ttp, it) also agrees 
with the DMRG. One additional prediction from the theory is 
that the pairing and density power laws have exponents that 
are inverse of each other While the DMRG estimates for the 
present system sizes do not satisfy this exactly, there are likely 
strong finite size effects, and we would like to revisit this with 
larger systems. 



APPENDIX B: LIMITING CASES IN THE J-K MODEL 

1. DBL[2,1] as Jordan-Wigner in a Model with No-Double 
Occupancy of Rungs 

In the DBL construction Eq. (|2]i, we can view one deter- 
minant as affecting a generalized flux attachment or Jordan- 
Wigner (JW) transformation and view the other determinant 
as a ground state of the new JW fermions. In the DBL[2,1] 
case, the d2 determinant is composed entirely from the bond- 
ing orbitals and becomes Eq. ( fT9] i; in particular, it prevents 
two bosons from being on the same rung. If we consider 
a hard-core boson model that prohibits double occupancy of 
rungs, then sign(dct2) is the conventional ID chain JW trans- 
formation on this ladder. As can be seen from the bottom 
panel in Fig.]?] our ring model in the DBL[2,1] phase appears, 
by its own dynamics, to strongly suppress double rung occu- 
pancy. The optimized power of dct2 in the trial wavefunction 
is also relatively small, making the sign(dct2)|dct2|^ factor 
look more like the Jordan-Wigner. Finally, the di determi- 
nant is composed of the bonding and antibonding orbitals, and 
its optimized power is relatively large though smaller than 1, 
which suggests that the JW fermions are not far from being 
free despite the large ring exchanges and the restricted rung 
occupancy. 

On the other hand, the DBL[2,2] would have some dou- 
ble occupancy of rungs and would be an example of a non- 
trivial JW, but unfortunately this phase appears to be unstable 
as described in Sec. lA 31 The instability happens because of 
special kinematic conditions satisfied here, which allow di- 
rect Josephson coupling and locking between enhanced bo- 
son modes. However, we expect that on iV-leg ladders with 

> 3, DBL phases will exist which cannot be described by 
a conventional Jordan-Wigner approach. 



2. Solution of the K-Only Model 

The pure ring model on the 2-leg ladder can be solved 
exactly. In the absence of boson hopping terms (J = 0), 
the number of bosons on each rung, A^rung(2^) = n{x, 1) + 
n{x, 2), is separately conserved. The number of bosons in 
each chain, iVch.i or -/Vch.ii, is also conserved. Any rung with 
^rung = or 2 effectively breaks the system into decoupled 
pieces. 

Consider an isolated segment of L rungs with each A^iung = 
1. The ring Hamiltonian is mapped to an XY spin chain by 
identifying boson configuration {n{x, l),n{x,2)} = {1,0} 
with spin up and {0,1} with spin down: 

L-l 

i7[L] = A'^[T+(a;)T-(a; + l) + H.c.] , (Bl) 

x = l 

where r are the usual spin- 1/2 operators. This is readily 
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solved by free fermions, and the ground-state energy is 



n=l 



sm 



-K 



7r(2rtinax+l) 
2(L+1) 



- 1 



(B2) 



(B3) 



2(L+1) 



where rimax = L/2 if L is even (in this case iVch.i = 
Nch.n = L/2), while 7i,nax = (i + l)/2 if L is odd 
(-/Vch.i = A^ch.ii ± !)■ In either case, the ground-state energy 
per boson is minimized by making L large, with the asymp- 
totic behavior 

_ E,.,[L] 2K Kil-2M 1 

epoiboson — r ~ ' ? ^ U(-pr {ti'l) 

L TT L L'^ 

Going back to the iC-only model on an infinitely long lad- 
der, we conclude that for arbitrary boson density p < 1/2, it 
is advantageous to phase-separate into an empty region and a 
half-filled region, since this minimizes the energy per parti- 
cle (p > 1/2 can be treated by particle-hole transformation). 
The phase separation arises because the ring exchanges pro- 
vide an effective attraction between particles - see also the 
next section. Note that the half-filled region itself is a highly 
correlated state of bosons. 



3. Bound State of Two Bosons for Large K 
and d-Wave Paired State at Low Densities 

Consider two bosons on an otherwise empty ladder In the 
absence of the hoppings, J = J± = 0, the states 



b{x, l)t6(x + 1, 2)t - b{x, 2)^b{x + 1, l)t 
V2 



I vacuum;,) (B5) 



are degenerate ground states with energy —K. Each can be 
viewed as a d-wave pair sitting on a plaquette. Small boson 
hopping can be included perturbatively, and the pair starts to 
hop along the ladder with amplitude J^/K. In fact, inspired 
by similar considerations in the ID t — J model in Ref. 
the zero momentum bound state can be written exactly for the 
general J — J± — K model, and is created by 



X n=l ^ ' 



fe(a;,l)t6(a; + n,2)1' (B6) 
5(a;,2)t&(x + n,l)^l . (B7) 



This is defined for K > 2J and has energy 

4J2 



-Abound — — A' — 



K 



(B8) 



the pair "size" is ^ = 1/ log[i^/(2J)]. 

Note that the bound state does not utilize any rung hopping. 
Also, it competes with a free state of the two bosons where 
each moves independently and where the total energy is 



assuming J, Jl > throughout. As an example, for Jj^ = 
the bound state has lower energy than the state with two 
free bosons when K > Kc[J± = 0] = 2 J, while Kc[J± = 
O.IJ] = 2.74 J and Kc increases with increasing J±. 

For a very small density of bosons, it is then natural to pro- 
pose that for K > Kc the bosons will form such tightly bound 
pairs which in turn form a Luttinger liquid. This is the picture 
of the d-wave paired phase. 

As discussed in the preceding section, in the absence of 
the hopping (or for very large K), the system at any den- 
sity phase-separates into empty and half-filled regions. This 
maximizes the ring energy per particle which attains —2K/ tt 
and, in particular, wins over the asymptotic — A'/2 energy per 
particle in the paired state. However, bosons in the phase- 
separated half-filled region form a non-trivial quantum state 
and can hardly gain any kinetic energy once the hopping is in- 
cluded. In the boson-paired liquid, on the other hand, the pairs 
can move and gain some kinetic energy, which can then stabi- 
lize the system against phase separation. As a crude estimate, 
we ignore the kinetic energy in the phase-separated case and 
compare — 2A'/7r with Eq. ( IB8I ) divided per boson and find 
that the boson-paired liquid is more stable for K < 3.8J. 
Thus, there is a sizable window over which the boson-paired 
liquid is energetically preferable over the free boson liquid 
and is stable against phase separation. 

At very low density we also need to consider more-boson 
bound states that successively improve their ring energy but 
lose kinetic energy, and the corresponding liquids could inter- 
vene between the boson-paired one and the phase separation. 
However, it seems likely that a sizable region of the boson- 
paired phase will remain; of course, it also competes with 
other phases such as DBL, as happens in Fig. |4]for the fill- 
ing p = 1 /9 studied by DMRG. 



APPENDIX C: TRIAL WAVEFUNCTIONS FOR 
BOSON-PAIRED PHASES 

Let us sketch how one can construct trial wavefunctions for 
boson-paired liquids such as the d-wave paired phase at low 
densities described in the last section and in Sec. lIIIDl and the 
s-wave paired phase described in Sec. IIIICI Let us denote the 
boson pair wavefunction as (f){r — r'). For a given pairing-up 
of particles we could write (t>{ri — r2)4>{r^ — Vi) . . . times 
some Jastrow-type factor for the center of mass coordinates 
Ri2 = (''I + '"2)/2,i?34 = ('"3 + ''4)/2, . . ., and then sym- 
metrize this to obtain a bosonic wavefunction; 



^ ^ S[(j){ri-r2)4>{r3-r4).. 
X Jastrow(i?i2, i?34, . . . )] . 



(CI) 
(C2) 



-2(2J + Jl 



(B9) 



Even if we take a simpler Jastrow which is already a symmet- 
ric function of {r^}, the symmetrization of the pairing part 
gives a permanent which is prohibitive in numerical calcula- 
tions. 

An alternative way to construct a symmetric bosonic wave- 
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function is 

* = Pi[g{n ~ r2)g{r3 - u) . . .] (C3) 
X det[ri,r2,r3, . . .] x Jastrow[ri, ra, rg, . . .] .(C4) 

Here pairing is realized with the help of a Pfaffian, which is 
a BCS wavefunction for spinless fermions and is specified by 
some pair function g{r—r'), while the total bosonic symmetry 
is recovered by multiplying by a second fermionic wavefunc- 
tion taken for simplicity to be a Slater determinant of some 
orbitals. Note that for tightly bound pairs when (j>{r — r') and 
g{r—r') are very short ranged the two constructions are essen- 
tially similar for statistically significant configurations. The 
Pfaffian and determinant are simple to work with in VMC, 
which allows one to construct trial wavefunctions also beyond 
the tightly bound limit. Of course, both the spinless fermion 
pair function g{r ~ r') m the Pfaffian and the orbitals in the 
determinant need to be judiciously chosen to realize the par- 
ticular boson pair function (p{r — r'). We should also note 
that such det x Pf wavefunction does not necessarily repre- 
sent a paired state of bosons, and one needs to measure the 
long-distance properties first. For example, we found that for 
a ID chain a product of a Slater determinant and a p^-wave 
BCS wavefunction has QLRO in the boson correlator when 
the -pairing problem is at weak coupling, while the boson 
correlator becomes short-range - as is needed for the boson- 
paired state - only when the p^; -pairing problem is at strong 
coupling. 

1. D-Wave Paired Wavefunction 

Consider for illustration low-density bosons in the Jj^ = 
limit. We can construct a d-wave paired state by taking 

9{xi,yi]X2,y2) = '5yi#y2sign(2/i - 2/2)e"''''"^"'/^ (C5) 

in the Pfaffian part, while combining the determinant and Jas- 
trow into one factor 

det X Jastrow = l\{ sin . (C6) 



The power of the determinant is understood as (det)P = 
sign(det)|det|^. The fermion pair function is p^-type and 
here requires that the two particles are on different chains 
(which is why it is strictly appropriate only in the = 
limit). For an isolated pair of bosons, this indeed gives the 
exact d-wave pair of Sec. IB 3 1 when p ^ 0. We use this wave- 
function with two variational parameters ^ and p to connect 
with the DMRG results at p = 1/9 in Sec. 1111 PI One could 
construct a similar boson-paired state by taking the Pfaffian 
factor to be a solution of a general spinless BCS problem with 
Py pairing on the ladder; this would allow more freedom (e.g., 
if we want to study Jj^ ^ 0), but has not been explored. 
2. S-Wave Paired Wavefunction 



We have also tried the det x Pf construction for the s- 
wave paired phase of Sec. 1111 CI Here we take the determinant 
to contain both bonding and antibonding orbitals, which re- 
sembles the deti in the DBL[2,1] or DBL[2,2] containing di 
fermions with dominant hopping along the chains (in partic- 
ular, the di bonding and antibonding bands become equally 
populated in the Jj^ = limit). We take the Pfaffian to be 
that for a p2:-wave spinless BCS problem on the ladder, which 
resembles d2 fermions trying to enter the DBL[2,2] regime 
but becoming paired. There are many parameters here such as 
the d2 intra- and inter-chain hoppings and px pairings. We 
have found that the boson s-paired phase from Fig. [3] and 
Sec. unci is reproduced roughly when the d2 BCS problem 
is in the weak-pairing regime for both the bonding and anti- 
bonding bands and has predominant inter-chain pairing. Note 
that the above di — d2 motivation of the det x Pf content is 
very loose and is hardly inspired by any gauge theory thinking 
but rather found by trial and error We also mention that there 
are likely other ways to construct such paired boson phases. 
For example, coming from 2D where paired bosons can form 
a true condensate, one natural construction to explore could 
be Pf X Pf . Nevertheless, the det x Pf that we studied gives 
us complementary confidence in the DMRG identification of 
the s-wave paired phase. 
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In this small system with 8 bosons, Fig.|7l there is also an abrupt 



change inside the DBL phase where the ground-state momentum 
jumps from (fc^, ky) — (0, 0) to {kx, ky) — {2tt/3, tt) at around 
K ~ 2.7 J. We can understand these quantum numbers from 
specific orbital occupations in the DBL wavefunction approach, 
and the change of the total ground-state momentum arises from 
a change of orbital occupations in the d2 bands. Near /\ ~ 2 
VMC finds the optimal DBL state to be deti [N[°^ = 8, iV^"' = 

0]abc X detaliV^"' = 6, A/'^"' = 2]abc- The label "abc" stands 
for antiperiodic boundary conditions when specifying the mean 
field orbitals; the total bosonic wavefunction satisfies periodic 
boundary conditions and carries momentum (0, 0). For the larger 
K, VMC finds the optimal DBL state deti[iVf' = 8, A^}"' = 

0] X det2[A^2°^ = 5, A^2'^' — 3], where now the orbitals are spec- 
ified with periodic boundary conditions; this state carries momen- 
tum (fcj;,fcj;) = (27r/3, tt). 

Note that AE in Fig. II If b) maintains a positive sign through- 
out the DBL phase, which can be understood roughly as fol- 
lows: A boson propagates via a bonding mode (with wavevector 
gi"') in one half and via an antibonding mode (qi^^) in the other 
half. The total phase accumulated by going around the system is 
(Qx'^ + qx'^)Lx/2 = TvpLx ~ mod 2it for p = 1/3 and the 
shown Lx multiples of 6. Thus, the zero twist boundary condi- 
tion provides the best (crude) matching for all K. When we tried 
different Lx (not shown), the AE could be of either sign. 



